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A search of axions produced in nuclear transitions was performed at the Kuo-Sheng Nuclear Power 
Station with a high-purity germanium detector of mass 1.06 kg at a distance of 28 m from the 2.9 GW 
reactor core. The expected experimental signatures were mono-energetic lines produced by their 
Primakoff or Compton conversions at the detector. Based on 459.0/96.3 days of Reactor ON/OFF 
data, no evidence of axion emissions were observed and constraints on the couplings g all and g a ee 
versus axion mass m a within the framework of invisible axion models were placed. The KSVZ 
and DFSZ models can be excluded for 10 4 eV < m a < 10 6 eV. Model-independent constraints on 
Qa^-glNN < 7-7 x 10" 9 GeV -2 for m a < 10 5 eV and g aea -gl NN < 1.3 x 10~ 10 for m a < 10 6 eV at 
90% confidence level were derived. This experimental approach provides a unique probe for axion 
mass at the keV— MeV range not accessible to the other techniques. 

PACS numbers: 14.80.Mz, 29.40.-n, 28.41.-i 



I. INTRODUCTION 

The axions (a) [l| were proposed in the 1970's as so- 
lution to the strong CP problem — the near-perfect can- 
cellations between the QCD vacuum angle and the quark 
mass matrix, as constrained by experimental bounds on 
the neutron electric dipole moments. The interactions of 
the axions with matter (photons, electrons and nucleons) 
can be described by the effective Lagrangian [l|, Q 



■^int — ga~)-f4>aE ■ B 
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where <p a , (E,B), ip e and yjpj 



represent respec- 



tively the axion, electromagnetic, electron and nucleon 
fields. The couplings of the axions to the photons and 
electrons are parametrized by ga-y-y and g a ee, while g^ NN 
and g^NN are their isoscalar and isovector couplings to 
the nucleons. 

A generic feature of the axion models is that all the 
coupling constants as well as the axion mass (m a ) are 
inversely proportional to the symmetry breaking scale 
(f a )- The original "Peccei-Quinn- Weinberg- Wilczek" 
(PQWW) model took f a to be the electroweak scale, 
implying m a of the order of ^100 keV. This has been 
tested and excluded after extensive efforts. Current re- 
search programs [4] focus on larger f a in the "invisible 
axion models" , the two popular variants of which are the 
DFSZ (or GUT) 0] and KSVZ (or hadronic) models @. 
The axion couplings with matter within the framework 



of these models were evaluated and discussed in details 
in Refs. 

The light-mass axion is a well-motivated dark mat- 
ter candidate. Cosmological and astrophysical ar- 
guments P, [t| constrain the axion mass to be 
10~ 6 eV < m a < 1CU 2 eV, but the bounds are model- 
dependent and with large uncertainties. Experiments 
have been performed to look for dark matter axions as 
well as those produced in the sun, power reactors and 
radioactive nuclear transitions. 

All previous reactor experiments EH focused on 
the searches of the PQWW axions via their decays, and 
contributed much to exclude their existence. In this ar- 
ticle, we present results on a new search over a broad 
axion mass range, using an alternative detection strat- 
egy through its interactions with matter. This detection 
scheme was successfully used in a previous experiment 
using radioactive isotope as axion source Q. 



II. REACTOR AS AXION SOURCE 

The axions are pseudoscalar particles and have quan- 
tum numbers like those of magnetic photons. It can 
be emitted through magnetic transitions in radioactive 



gamma-decays jl2J. Nuclear power reactors are power- 



ful radioactive sources and are therefore potential axion 
sources as well. Axions can be emitted in competition 
with the photons as a result of neutron captures 
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TABLE I: A summary of magnetic transitions and their esti- 
mated fluxes at a typical 2.9 GW power reactor. 
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The axion lifetime is then given by 

1 
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or nuclear de-excitations 
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There are six prominent channels of magnetic gamma- 
transitions at typical nuclear reactors, as listed in Tabled] 
Thermal neutron captures on the 10 B in the control rods 
and on proton in the cooling water produce a+ 7 Li* and 
d+7, respectively. Their photon fluxes ($7), in units of 
fission -1 and GCLwere evaluated by full neutron trans- 
port simulations The other sources of 91 Y*, 97 Nb*, 
135 Xe*, and 137 Ba* are all fission daughters. Their cor- 
responding I, were derived from standard tables on fis- 



sion yields 1J] and cross-checked by previous calcula- 
tions For comparisons, the v e 15] and i/e-yields 
at reactors are about 7.2 and ~ 10 -3 fission -1 , respec- 
tively. 

The axion flux (4> a ) at a distance L from a reactor core 
of fission rate Rf can be described by 



III. AXION DETECTION AND EXPERIMENTAL 
SETUP 

Data were taken with a high-purity germanium detec- 
tor (HPGe) of mass 1.06 kg at the Kuo-Sheng (KS) Re- 
actor Laboratory. The HPGe target and the associated 
anti-Compton (AC) detectors as well as passive shieldings 
are depicted in Figure [T] The principal AC detector was 
a well-shaped Nal(Tl) crystal scintillator of mass 19.7 kg. 
Descriptions of the experimental hardware and analysis 
procedures were presented in details in Refs. [igfejlq]. 
The primary scientific goal was the search of neutrino 
magnetic moments. A physics threshold of 12 keV 
and a background level of 1 day -1 kg _1 keV -1 compa- 
rable to those of underground dark matter experiments 
were achieved. The source-detector distance of KS was 
L=28 m while p^ was modeled to be 27.75 m of water 
and 0.25 m of lead. 



ML) = ^^■^■P dk -Pint , (4) 

where (r a /T 7 ) is the branching ratio of axion emissions 
in the transitions. It depends on the axion-nucleon cou- 
plings and the nuclear structures of the transitions. The 
probabilities of the axions surviving the flight from re- 
actor core to detector without decays or interactions are 
given, respectively, by 

P dk = exp[-:^] (5) 

Pa ■ T a 

and 

Pint = exp[-L ■ p L ■ <r int ] , (6) 

where m a , T a , p a , E a are the axion mass, lifetime, mo- 
mentum and total energy, ai n t is the axion interaction 
cross section with matter at effective target number den- 
sity p L . 

In particular, the axions can decay in flight via the 
emissions of 27 (r 77 : a — > 77) or e + e~ pairs (r ee : a — > 
e + e - ). Their decay rates are related to the ga-yy and g aee 




FIG. 1: Schematic layout of the HPGe with its anti-Compton 
detectors as well as inner shieldings and radon purge system. 

The search strategies for reactor axions with these 
unique HPGe data were inspired by a previous experi- 
ment Q where a 15 kCi 7-source of 65 Zn was used as 
a potential axion source instead. Two interaction mech- 
anisms of axions with matter were studied: Primakoff 
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FIG. 2: Schematic diagrams of interactions of axions with 
matter, via (a) Primakoff and (b) Compton conversions. 



The total energy of the axions can be fully converted 
by either of the two processes into measurable ioniza- 
tion energy in the HPGe, such that the experimental sig- 
natures are the presence of mono-energetic lines at the 
known E a 's during the Reactor ON periods. In compar- 
ison, previous reactor-based axion experiments studied 
instead the axion decay channels T 77 (lp| or T ee 11|. 
They were therefore not sensitive to the invisible axion 
regime where m a are very small and decays are kinemat- 
ically blocked or suppressed. 



and Compton conversions as shown schematically in Fig- 
ures [2h&b. These processes are respectively and inde- 
pendently sensitive to g a ~n and g aee . Their cross-sections 
were both given in Ref. [2J — (a) Primakoff conversion on 
the nuclei: 



IV. DATA ANALYSIS 
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where a em is the electromagnetic coupling, Z is the 
atomic number of the target, f3 — p a /E a and x is the 
atomic-screening correction factor given in Eq. 20 of 
Ref. 0; (b) Compton conversion on the electrons: 
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where y = 2m e E a 



The m a dependence of the two 



cross-sections at the pn — > d7 transition energy of E a = 
2.23 MeV, using the normalizations of g a77 = 1 GeV _] 
and g aee = 1, arc illustrated in Figure [3] 
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FIG. 3: The variations of the Primakoff and Compton con- 
version cross-sections with the axion mass at E a = 2.23 MeV, 
using the normalizations of g all = 1 GeV -1 and g aee = 1. 
The dashed line represents the case where the atomic screen- 
ing effects are switched off (x — !)• 



The signal rates (Sp/c) for axion Primakoff/Compton 
Conversions in one kilogram of target mass are given by 



S P = op ■ 4> a ■ N ■ ep 
and 

Sc = o c ■ (j) a - Z ■ N ■ e c 



(12) 
(13) 



where e P / c are the efficiencies of full energy deposition 
at the HPGc detector, N is the number of atoms in the 
kilogram target, and Z accounts for the electron target 
number in Compton process. The various efficiency fac- 
tors were evaluated by full simulations and listed in Ta- 
ble [TTJ Full energy depositions for Compton conversion 
at the HPGe detectors were due to interactions only in 
Ge, such that only iV(Ge), Z(Ge)=32 and ec(Ge) were 
involved in the derivation of Sc- However, photons from 
Primakoff conversions in both Ge and Nal could con- 
tribute to Sp as full-energy peaks at the HPGe, such that 
there are terms involving, respectively, [7V(Ge), ep(Ge)] 
and [iV(Nal), e P (NaI)]. 

Evidence of reactor axions would manifest as peaks at 
the known energies of Table U in the Reactor ON— OFF 
residual spectra in HPGe. Following the naming con- 
ventions of Ref. Periods-I (July 2001- April 2002) 
and -III (Sept. 2004-Oct. 2005) with 180.1/52.7 and 
278.9/43.6 days of the KS-HPGc Reactor ON/OFF data, 
respectively, were used for analysis. Candidate events 
were those uncorrelated with the anti-Compton and 
Cosmic-Ray vetos and having pulse shapes consistent 
with 7-events. Selections of these events and their ef- 
ficiencies were discussed in Ref. [lj|. As illustrations, 
the Period-Ill ON/OFF background spectra and the 
ON— OFF residual spectra for the six candidate lines are 
depicted in Figures 0h&b, respectively. The background 
7-lincs were identified 15j and indicated that ambient 
radioactivity dominated. 

The count rates and their errors of the various transi- 
tions shown in Table [TT| were derived by best-fits of the 
residual spectra to Gaussian lines at fixed E^s and res- 
olutions. No excess were observed in all channels and 
upper limits of the signal rates (S u ) at 90% confidence 
level(CL) were derived. The most sensitive channel in 
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TABLE II: Summary of the detector efficiencies for Primakoff and Compton conversions on axions of energy E a , the measured 
Reactor ON— OFF event rates at the signal regions and the 90% CL upper limits (S u ). 
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FIG. 4: (a) The after-cut ON and OFF spectra and (b) the 
residual spectra for the six specific channels for Period-Ill, 
with the best-fit Gaussian peaks overlaid. 



terms of S u is the 2.23 MeV transition (RMS resolution 
2.3 keV) in the np^da interaction, because of the lower 
background level compared to those at ^500 keV (by ~ 
10" 2 [H). 

For completeness, searches were also performed 
at: (a) half E a to look for single-7 absorption in T 77 de- 



cays, and (b) full E a of the individual ON/OFF spectra 
to look for axion emissions from other possible steady- 
state sources, such as those from the sun [13, EH- No 
signals were observed in both cases. 



V. DERIVATIONS OF AXION PARAMETERS 

The experimentally measured upper limits S u of Ta- 
ble pi] can be translated to bounds among the axion pa- 
rameters: m a , g Q77 , g aee , g° NN and g\ NN - Two ap- 
proaches were adopted: (a) specific models were used and 
tested, and (b) model-independent constraints among the 
parameters were derived. 



A. Branching Ratios for Axion Emissions 

The limits S u are related to reactor axion emission and 
detection via 



Sp 



(14) 



Both Sp and Sc depend on the reactor axion flux, and 
thus the branching ratios (r a /T 7 ). 

Starting from the interaction Lagrangian of Eq. [T] the 
axion emission branching ratio of the pn — * d7 isovector 
Ml transition can be expressed as T3, 19j : 
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In the formulae, L is the multi-polarity of the transition, 
5 is the known E(L + 1)/ML mixing ratio, /iq and 



5 



are respectively the isoscalar and isovector magnetic mo- 
ments which can be derived from the proton and neutron 
magnetic moments (fi p and /j. n ) via: 



Mo = M P + A*n — 0- 88 
and 

Mi = M P - A*n - 4.71 



(17) 



(18) 



The nuclear physics of the transitions are parametrized 
by the matrix elements (3 and r\ defined as 



TABLE III: The calculated numerical factors A 0/1 of Eq. [25] 
under the KSVZ and DFSZ invisible axion models using the 
formulae of Ref. [7(. 
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A 
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(20) 



where Ji and Jf are the initial and final nuclear angular 
momentum in the transitions, while and a(i) are the 
orbital angular momentum and nuclear spin operators. 



B. Invisible Axion Models 

The evaluations of (r a /r 7 ) given in Eq. [TB] involve 
modeling of the couplings g^ NN and g\ NN . Within the 
framework of the invisible axion models, these couplings 
are inversely proportional to the symmetry breaking scale 
f a which in turn is related to m a 0, Q via 



m a /eV = ( ,,„,. )[- 



fa/GeV 'A + z' 



(21) 



such that the couplings depend linearly on m a . 

The formulae given in Table 1 of Ref. were adopted 
for the parametrizations of these axion-nucleon cou- 
plings: 



n° - C 

9aNN — ° 



(3F-D) 
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{X u -X d - N f ) 



+ -(X u + 2X d - Nf) 



and 



1 - z„ 



9aNN ~ ° ' 2 ^ " _ d | l. - 

where the factor common to both couplings is 
C = 5.2 x 1(T 8 (JL)[ mo /eV] . 

IV f 



(22) 



(23) 



(24) 



The terms D ~ 0.77 and F ~ 0.48 [l8| are the reduced 
matrix elements for the octet axial vector currents, S — 
0.33 ± 0.04 [20| denotes the flavor singlet axial charge, 
Nf = 3 is the number of families, z = (m u jm,d) — 0.56 
is the ratio of the up-to-down quark mass, while X u and 
Xd represent respectively the PQ charge of the u and d 
quarks. The KSVZ and DFSZ models differ essentially 
in their choices of (X u , Xd). 



Previous evaluations of axion fluxes in the stellar [21 1 
and solar 17 , IH environment adopted the KSVZ model 
where X u = Xd = 0. This model also specifies g ae e = 
at tree level such that the results are only applicable to 
probe the axion-photon g a77 -couplings. We extended the 
analysis to include also the DFSZ model, which allows 
finite g a77 - and g Qee -couplings. The parameters X u and 
Xd are positive-definite constrained by X u +X d = 1. The 
values of X u = Xd = 0.5 were chosen for this analysis. 
Defining 



o/i 

9aNN 



A ' 1 [m Q /eV] 



(25) 



the calculable numerical factors A / 1 under both KSVZ 
and DFSZ models are tabulated in Table Hill The isovec- 



tor couplings g\ NN are the same for both models with 
this specific choice of (X u , Xd). 

Once g® NN and g aNN are fixed by the invisible axion 
models, the evaluations of (r a /r 7 ) depend on the nu- 
clear physics inputs: 5, (3 and r\. The values adopted for 
these parameters are summarized in Table IIVI Among 
the transitions, the pn — > d7 is a pure isovector Ml pro- 
cess. Its branching ratio (r a /r 7 ) p „, as given by Eq. [TBI 
is independent of (S, /3, 77) and was used in previous reac- 
tor axion experiments [13, EH, [H| • The 7 Li* transition 
is also pre-dominantly Ml. The neutron shell is closed 
and the transition is driven by the odd-proton, such that 



(S,P,rj) ~ (0,1,0.5) was adopted as in Ref. [17[. There 
were no calculations on (/?, 77) for the remaining four fis- 
sion daughter isotopes. To make estimations, we project 
from the results on other heavy isotopes [2j, [21(, also 
summarized in Table IIVI The matrix elements (j3, rf) 
for heavy isotopes with unpaired proton and neutron are 
taken to be ~ (1, —3) and ~ (—1, 1), respectively. 

Once these (/3, rf) assignments are made, the branching 
ratios (r a /r 7 ) can be readily evaluated with Eqs.[T5Tl2~51 
In particular, j3 < for the odd-neutron nuclei 135 Xe* 
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TABLE IV: A summary of the nuclear physics input in the evaluations of (r a /r 7 ). Calculated values of (/3, rj) on the other 
isotopes from previous works are included for comparisons. The branching ratios (F a /F 7 ) at m a ~l eV under the KSVZ and 
DFSZ invisible axion models and the quality factors (QF) relative to that for the pn — * d"/ channel are also shown. 



Channel Unpaired Transition 

p/n 



/3 



Remarks 



(r„/r 7 ) t QF/(QF) pn 

KSVZ Model 



-TTT 

-15 
-15 
-15 
-16 

-TT- 

-16 
-16 
-15 



(r a /r 7 ) t qf/(qf) p „ 

DFSZ Model 



-nr 

-15 
-15 
-15 
-18 

-18 



pn — > d7 
7 Li* 

91 y* 

97 Nb* 

135 Xe , 



137 



Ba* 



57 Fe* 
55 Mn* 
Na* 



P 
P 
P 
n 
n 



Ml 
Ml 
M4 
M4 
M4 
M4 






1 


0.5 





1 


-3 





1 


-3 





-1 


1 





-1 


1 



Refs. [10, 11,12} 
Ref. [17] " 

Inferred 
{ from 
Refs. [2, 21] 



7.4 x 10 

4.0 x 10 

2.1 x 10 
2.1 x 10 
1.6 x 10 
1.6 x 10 



1.00 
0.50 
0.33 
0.13 
0.02 
0.01 



7.4 x 

3.0 x 

1.5 x 
1.5 x 

7.1 x 
7.1 x 



1.00 
0.43 
0.28 
0.11 
0.003 
0.003 
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Ml 
Ml 
Ml 
Ml 



0.44 1.81 -6.59 Ref. [2] 

0.002 -1.19 0.80 Ref. [21] 

0.052 0.79 -3.74 Ref. [21] 

0.058 0.88 -1.20 Ref. [21] 



1.0 x 10 

3.4 x 10 

8.5 x 10 
1.9 x 10 



7.1 x 

2.2 x 
6.6 x 
1.4 x 



1Q -io 

io- 17 

1Q -16 

IO" 15 



Evaluated at m a =l eV 



and 137 Ba*, such that the g^ NN and g\ NN terms have op- 
posite signs and (r a /r 7 ) can be vanishingly small within 
the large nuclear physics uncertainties. These two chan- 
nels were discarded in subsequent analysis. The varia- 
tions of (r a /T 7 ) with m a of the four remaining channels 
are depicted in Figure[5l The differences within the chan- 
nels and between the KSVZ and DFSZ models are small 
relative to the scale in the log-log plot, and are repre- 
sented by the width of the line. 
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log, [ m.(eV) ] 

FIG. 5: The variations of (r a /T 7 ) with m a in both KSVZ 
and DFSZ invisible axion models as parametrized in Ref. Q. 
The 135 Xe* and 137 Ba* channels have large uncertainties and 
were discarded. The width represents the variations within 
the four remaining channels and among the two models. 

The experimental sensitivities in g aT y and g aee can be 
described by a quality factor (QF) which is related to the 
7- flux $ 7 of Table fl] and the upper limit signal rates S u 
of Table 02 via: 

^ aQFs S]p5. (26) 

Both (r a /r 7 ) and the relative QF (with respect to that of 



the pn — > d7 channel) evaluated at m a =l eV were shown 
in Table IIVI to illustrate the relative strength of the vari- 
ous channels. It can be seen that the leading contribution 
to the sensitivities is from the pn — > d7 channel. 

C. Model-Dependent Limits 

The (r a /r 7 ) estimates in Tabic [TV] for pn — > d7 and 
7 Li* are accurate while those for the four heavy fission 
isotopes are expected to have large uncertainties. In ad- 
dition, the relative QF values of Table IIVI indicate that 
the pn — > d'j channel mostly defines the sensitivities. Ac- 
cordingly, we take the conservative approach that the ex- 
perimental limits on <7 a77 and g aee were derived only from 
this channel. The two detection channels were treated in- 
dependently, via the relations 

S P <S U & S C <S U , (27) 

assuming, respectively, g aee = and g ail — 0. These 
produce less stringent bounds compared to those from 
the convoluted case of Eq. [Jj] such that the results have 
general validity. 

The exclusion plot of g ail versus m a is depicted in Fig- 
ure [5] The model-dependence was introduced by fixing 
(r a /r 7 ) p „ as a function of m a via Eqs. [T5&l2"5l The val- 
ues of (T a /Tj) pn depend only on g\ NN and are therefore 
the same for both KSVZ and DFSZ models, as shown in 
Table HVl Consequently, both models produce a common 
exclusion region, denoted by "KS". On the other hand, 
the KSVZ model specifies g aee — at tree level, such 
that only the DFSZ model can be meaningfully applied 
to define the exclusion region, also labeled as "KS", in 
the g ae e versus m a plot of Figure [7] 

In both Figures [SHZl the vertical bounds labeled 
"Kine" at 2.23 MeV are due to the kinematical con- 
straints from the maximum E a . The sensitivities are sup- 
pressed at "Decay" for the large (m a ,<7 a77 / 1 g a ee) regions, 
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FIG. 6: Exclusion plots of <? a77 versus m a for g a ee~0. The 
limits from the KS experiment at 90% CL are denoted by 
"KS". They are derived by fixing (r a /F 7 ) pn with the KSVZ 
and DFSZ invisible axion models. Predicted regions of the 
PQWW, KSVZ, DFSZ and HW models on the (m a ,c/ a77 ) 
plane are overlaid. The boundaries denned by "Int", "De- 
cay" and "Kine" are constraints due to axions interactions 
with matter, decays in flights and kinematics, respectively. 
The bounds marked "Zn" are results from Ref. ||, while re- 
gion labeled "R" are from previous reactor axion experiments 
studying T 77 Results from the other axion experiments 
using different techniques [23j, 24, 25, 26, 29] are displayed as 
colored blocks. The astrophysical bounds 9] are denoted by 
the striped region. 



due to r 77 and T ee decays in flight. The lack of sensitiv- 
ities at "Int" for large g a ~ n > 20 GeV -1 and g aee > 0.2 
are due to axion interactions in the matter between reac- 
tor core and detector. Limits marked "R" and "Zn" are 
respectively from previous reactor experiments studying 
r 77 [10j or r ee ll| and the radioactive source experi- 
ment [2j. They were derived using the same modeling 
schemes on (r a /r 7 ) as "KS". The bounds from the KS 
reactor axion searches improve on those of Ref. [2| by two 
orders of magnitude, owing to enhanced axion flux, lower 
background and larger data sample. 

The KS results define the global exclusion boundaries 
in g ail for 10 3 eV < m a < 10 6 eV and in g aee for 
10 2 eV < m a < 10 6 eV. Astrophysics arguments on 
stellar cooling and red giant yields 0| provide more 
stringent bounds for m a < 10 4 eV but these are model- 



FIG. 7: Exclusion plots of g aee versus m a for <7 a77 =0. Sim- 
ilar conventions as Figure [6] are adopted. The "KS" bounds 
are derived by fixing (r a /r 7 ) pn with the DFSZ model. Re- 
gions marked "R" was excluded by experiment studying r ee 
at reactors [TlJ. Bounds from othe experiments [23L 1271 . |2fi| 
are shown as colored blocks. 



dependent. They are represented by the striped regions 
in Figures respectively. Comparisons of the KS ex- 
cluded regions with the KSVZ/DFSZ predictions on the 
(m Q , ga^j/gaee) planes would rule out these models at 
10 4 eV < m a < 10 6 eV. An example of other existing 
model pre dicting axion mass at the MeV range is the HW 
model [22|, also depicted in the figures. 

The experimental approach presented in this article 
can probe the keV— MeV axion mass range which is not 
accessible to the other techniques. At large m a > 10 6 eV, 
the sensitivities in both g a ^ and g ae e are dominated 
by the accelerator-based "beam dump" experiments [23| • 
Exclusion boundaries at small m a are defined by: (a) 
for g ail — the germanium [24| and CAST [25[ experi- 
ments studying solar axions, and the axion dark matter 
searches with microwave cavity [26j |; and (b) for g aee — 
the positronium decay [ilj and macroscopic force (28j 
experiments. At g ajl > 10~ 2 GeV -1 , the solar axion 
experiments are limited by axion interactions inside the 
Sun Part of this large-g a 77 region has been rejected 
by the optical laser experiments [29(, while the KS re- 
sults contribute to probe and exclude a remaining hole 
at m a ~ 10 1 eV. For completeness, we mention also the 
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recent PVLAS experiment which reported a finite light 
polarization rotation in vacuum with a transverse mag- 
netic field [3(|. This result was interpreted as the region 
"P" of finite (m a ,g a71 ) in Figure well contradicted by 
many other experiments. This would imply that it may 
not be appropriate to analyze the PVLAS results using 
the existing axion models. 



D. Model-Independent Constraints 

The KS exclusion regions of Figures [5$43 were evalu- 
ated within the framework of the invisible axion mod- 
els. Alternatively, it is also instructive to derive model- 
independent constraints among the axion parameters. 
Following the reasonings of Section V.C, only the results 
from the pn — > dj channel were adopted, and the two 
couplings g aTy and g ae e were treated independently. 

The experimental sensitivities are defined by the sig- 
nal rates S P / C which are proportional to 5a 77 / aee (r a /r 7 ) 
and alternatively, via Eq.[T51 to g 2 all / aye gl N N ■ Limit on 
the Primakoff conversion rate Sp < »!? u (pn — > d7) gives 
rise to the model-independent constraint 



(a) 



r 3a 77 ' ( r° )pn 

^ 1 

5a77 ' 9aNN 



< 5.9 x 1CT 17 GeV~ 2 (28) 

< 7.7 x 1(T 9 GcV" 1 1 (29) 



which are applicable within the ranges of m a < 10 eV 
and 3,277 < 20 GeV^ 1 , bounded by kinematics and axion 
interactions in flight, respectively. Similarly, the limit on 
the Compton conversion rate Sc < S u (pn — » d7) leads 
to 



{^) pn < 1.7x10- 



Jaee ' 9 a NN 



< 1.3x10 



-10 



(30) 
(31) 



for m a < 10 6 eV and g ae e < 0.2. The loss of sensitiv- 
ities in (7077 at m a ~ 10 6 eV can be explained by the 
reduction of the Primakoff cross-section as depicted in 
Figure [3] These constraints on <? a77 and g aee as func- 
tions of g\ NN and (r a /r 7 ) p „ are illustrated respectively 
in Figures. EK&b. Limits from the most sensitive labora- 
tory experiments are also displayed for comparisons. 

The KS limits on g a77 are not as sensitive as those of 
CAST [25| and the solar-germanium [24| experiments for 
m a < 10 4 eV at all branching ratios. On the contrary, 
the q a pj, sensitivities exceed those of the positronium de- 
cay [27j and macroscopic force 28[ experiments. New 
regions are probed for m a > 10~ 6 eV with (r a /T 7 ) p „ > 
10~ 9 and for m a < 10~ 6 eV with (r a /T 7 ) p „ > 10" 5 . 



VI. SUMMARY AND PROSPECTS 

This article reports the first study of possible emis- 
sions of axions from power reactors using Primakoff and 
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FIG. 8: The model-independent exclusion regions of the KS 
experiment for (a) g all and (b) g aee on the (r a /T 7 ) pn and 
9aNN axes. The g all and g a ee limits from the leading labora- 
tory experiments are also shown for comparisons. The ranges 
of validity in m a are indicated. 



Compton conversions as the detection mechanisms. No 
evidence were observed and constraints on axion param- 
eters were placed. 

The exclusion regions in the (m a , g a -yy/ gaee) param- 
eter space were identified within the framework of the 
invisible axion models. The branching ratios for ax- 
ion emissions associated with radioactive 7-decays are 
proportional to m 2 ,, such that the experiment is sensi- 
tive mostly at the large axion mass region. The KS 
results define the global exclusion boundaries for g a ~ n 
and g a ee and excluded the KSVZ and DFSZ models for 
10 4 eV < m a < 10 6 eV. 

Independent of models, the KS results are not as sensi- 
tive in constraining (^77 compared to those from solar ax- 



9 



ion searches, but improve on the limits in g aee for branch- 
ing ratios as small as (r a /r 7 ) p „ > 10 -9 for light-mass 
axions such as those within the cosmologically-preferred 
range of 10~ 6 eV < m a < 10~ 2 eV. 

Our studies therefore indicate that this experimental 
approach can provide competitive sensitivities compared 
to the other techniques when (a) the axion physics is 
correctly described by the invisible axion models and 
the axion mass is at the keV— MeV range, or (b) the 
axion physics allows relatively large axion-nucleon cou- 
plings (and consequently the axion emission branching 
ratios) at small axion mass not yet covered by the cur- 
rent theoretical modeling. 
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